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Abstract 

The Lamb shift results from the couphng of an atom to vacuum fluctuations of quantum fields, 
so corrections are expected to arise when the spacetime is curved since the vacuum fluctuations are 
modified by the presence of spacetime curvature. Here, we calculate the curvature-induced correc- 
tion to the Lamb shift outside a spherically symmetric object and demonstrate that this correction 
can be remarkably significant outside a compact massive astrophysical body. For instance, for a 
neutron star or a stellar mass black hole, the correction is ~ 25% at a radial distance of AGM/c^, 
~ 16% at lOGM/c^ and as large as ~ 1.6% even at lOOGM/c^ , where M is the mass of the object, 
G the Newtonian constant, and c the speed of light. In principle, we can look at the spectra from 
a distant compact super-massive body to find such corrections. Therefore, our results suggest a 
possible way of detecting fundamental quantum effects in astronomical observations. 
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The interplay of quantum theory and general relativity has brought new insights into our 
understanding of these two fundamental theories which are pillars of modern physics. The 
Hawking radiation of black holes is one of the most striking phenomena that are revealed by 
our attempt to establish quantum field theory in curved spacetime. On one hand, as one of 
the quantum effects that are unique to curved spacetime, the Hawking radiation, which has 
been extensively examined in many different contexts, has been considered as a " Rosetta 
stone" to relate quantum theory, general relativity and thermodynamics and is expected to 
be a indispensable part of a yet-to-be-found full theory of quantum gravity; on the other 
hand, it still remains a tremendous task to experimentally verify the Hawking effect as well 
as other quantum effects unique to curved spacetime, although a lot of efforts have been 
made to observe it in analogue systems [l| . In consideration of the lack of direct experiment 
test of the quantum effects unique to curved spacetime, here we take a different approach, 
i.e., we ask, what are the corrections, caused by the fact that the spacetime is curved, to 
quantum effects already existing in flat spacetime? and can these corrections be observable? 

In this regard, the quantum effect that first comes to our mind is the Lamb shift, a subtle 
difference in energy between two levels of an atom which arises as a result of the coupling 
of the atom to vacuum fluctuations of quantum fields. The Lamb shift is one of the most 
important and inspiring discoveries of the last century that marked the beginning of rnodern 
quantum electromagnetic field theory. Since its discovery in experiment in 1947 [2], the 
Lamb shift has attracted a great deal of attention and been measured with great precision. 
It is worthwhile to note that the measurement of the Lamb shift also plays an important role 
in our understanding of nuclear structure and in the determination of fundamental physical 
constants jsj-Q]. So far, on the theoretical front, the Lamb shift has been examined in various 
circumstances, for example, in the presence of cavities [8|, or in a thermal bath B-lli. It 



has recently been shown that the non-inertial motion of the atom also induces corrections 



to the Lamb shift [13l4l6l|. However, all the aforementioned studies are concerned with fiat 
spacetimes. Therefore, it remains interesting to see what happens if the atom is placed in 
a curved spacetime rather than a flat one. Now, a correction to the Lamb shift as opposed 
to its original value in a flat spacetime is generally expected, since vacuum fluctuations of 
quantum fields are modified by the presence of spacetime curvature ^. Here, we are concerned 
with the correction to the Lamb shift which is caused by the modified vacuum fluctuations 
due to the spacetime curvature outside a spherically symmetric astrophysical object, and, 
surprisingly, we find that this correction is potentially observable, thus providing a possible 
way of checking fundamental quantum effects in astronomical observations. 

The Lamb shift has been studied in different physical contexts^U yielding the same result 



in agreement with experiment with remarkable precision 13|, ll7H2l| . Our calculation here is 



based upon the so-called Dalibard, Dupont-Roc, and Cohen- Tannoudji(DDC) formalism 22 



231] in which the contribution of vacuum fluctuations to the energy-level shift of an atom 



^ Let us also note that the energy-level structure of an atom will in general be different from that in a flat 
space since it is now determined by the wave equation written in a curved space. 
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and that of the radiation reaction are distinctively separated. Since we are interested in the 
correction induced by the curvature, we consider, for simphcity, a two-level atom interacting 
with quantized real massless scalar fields in vacuum at a fixed radial distance outside a 
spherical massive body. When a curved spacetime is concerned as opposed to a fiat one, a 
delicate issue then arises as to how the vacuum state of the quantum fields is specified. In 



this paper, we assume that the scalar field is in the Boulware vacuum [2J], since it is the 
vacuum state outside a massive spherical body which has not collapsed through its event 
horizon. Two stationary states of the atom are represented by |+) and |— ) respectively. The 
Hamiltonian H of the system (atom+ bath of fluctuating scalar fields) is composed of the 
following three parts 

Ha{t) = hujoS-Xr) , (1) 

/dt 
d^khuj^ala^,-, (2) 

Hiir) = /i^2(r)0(x(r)) . (3) 

Here Ha{t), Hpij) and Hj{t) are respectively the Hamiltonian of the atom, the field and 
the Hamiltonian describing the interaction between them. 5*2 (0) = |^(|— )(+| — |+)(— I), 
'S'2(0) = |(|+) (+1 — |— ) (— I), and (p is the field operator, at and are the field creation and 
annihilation operators, fi is the coupling constant which is assumed to be small. 

Assuming the initial state of the atom is \b) and the field is in the Boulware vacuum, one 



can show, in the framework of Dalibard, Dupont-Roc, and Cohen- Tannoudji (DDC) 13 



22 



23 



that the contributions of vacuum fluctuations and radiation reaction to the energy shift 



of level b are respectively given by, 

,2 /-T 



{SE,)rr = -4- r dr'x^{x{r),x{r')){C\{r,r') , (5) 



TO 



where r is the proper time of the atom and x{t) represents its stationary trajectory which we 
will specify later. and are the symmetric correlation function and linear susceptibility 
of the field defined as 

C^(x(r),x(r')) = ^(O|{0^(x(r)),0^(x(r'))}|O) , (6) 

X^(x(r),x(r')) = i(O|[0/(x(r)),0/(x(r'))]|O) , (7) 
and {C'^)i, and {X^)h are those of the atom, which are given by 

{C\{r,r') = i^|(6|52(0)M)r(e^--^^ + e-'--^-) , (8) 

d 

{X\{r.r') = i5^|(6|52(0)|d)p(e--^--e— ^0 ■ (9) 



d 
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Here um = u^b — <^d is the energy spacing between the two levels, \ubd\ = <^o for b d, 
At = t — t' and the sum over d extends over a complete set of states of the atom. 

Now let us calculate the energy-level shift of a two-level atom at a fixed radial distance, 
of which the trajectory is described by 

1 



t(r) 



(r-ro), r(r)=r, e{r) = e, 0(r) = . 



(10) 



First, we need to solve the Klein-Gordon equation in the given curved background, define 
the Boulware vacuum and calculate the statistical functions of the field. For this purpose, 
let us note that the metric of the spacetime outside a spherical object can be written in the 
Schwarzschild coordinates as 

-1 



ds' 



^-^]dt' 
r 



dr' - r' {de' + sin' e d^') . 



(11) 



Here is the Schwarzschild radius, = 2GM/c^ in which M is the mass of the object, G 
is the gravitational constant and c the speed of light. Solving the Klein-Gordon equation 
satisfied by a massless scalar field 25|| , one finds, in the exterior region of the massive body, 
a complete set of normalized basis functions 



h 



Atiuc 
Ti 



Ri{uj\r) Yim{9,ip) 



(12) 
(13) 



where Yim{0, f) are the spherical harmonics and the radial functions satisfy the following 
differential equation 



d 



dr 



d 



dr 



c^(l — rg/r) 



1(1 + 1) 



Ri{oo\r) 







This equation can also be written as 

d^ 



dr* 



rRi{uj\r)) + [uj^ - Vi{r)]{rRi{Lo\r)) = 



with 



V,(r) 



1 



2M 



/(/ + !) ^ 2M 



(14) 



(15) 



(16) 



being the effective potential. Eq. (]T5l) admits the following asymptotic solutions, since 
VJ(r) ~ in the asymptotic regions (r ~ 2M and r — )■ oo) 



Riiulr) ~ 



Ri{u\r) 



r~^e'^'*l^ + Ai{u})r-^e-''^'''l\ r ^ r, , 
5/(w)r-ie^^''*/^ r ->cx) , 

fi;(a;)r-ie-^'^^-/'=, r ^ r, , 



(17) 



(18) 
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Here r^, = r + rsln(;^ — 1) is the Regge- Wheeler tortoise coordinate, and Bi^u) and Ai{u) 
are respectively the transmission and reflection coefficients, which obey 



\Muj)\ = \Muj)\, (19) 

i-\Mu)\' = i-\M^)\' = mco)\' . (20) 

The physical interpretation of these modes is that u represents the modes emerging from 
the past horizon and u denotes those coming in from inflnity. The Boulware vacuum state 
is deflned by requiring normal modes to be positive frequency with respect to the Killing 
vector d/dt with respect to which the exterior region is static 24 . 

The fleld operator can be expanded in terms of the complete set of the basis modes, 
and the Feynamn propagator (for t > t') of the massless scalar fleld in the vacuum is given 
by HE 



It then readily follows that C^{x{t),x{t')) = |(D^(x, x') + D^(x', x)) and x^{x{t),x{t')) = 
^{D^{x,x') — D^{x',x)). Inserting these statistical functions into Eqs. ([4]) and (E]), we can 
separately calculate the contributions of vacuum fluctuations and radiation reaction to the 
shift of an energy level of the atom. Adding them up, we obtain the total shift of level b 

2 /*oo 

SEt = -j^Y.\(b\S,md)n x[g(a;|r) + g(c.|r)], (22) 



d " V900 



where 



^ oo 

g{u:\r) = + l)\RMr)\\ (23) 



goo = l-rjr and EJ(&I^2(0)M)P = 1/4. 

The Lamb shift of the atom is obtained by subtracting the shift of the ground state from 
that of the excited state 

^ r dto 

^ ~ ~64^ 

with 



^'^ ^ -/goo " Vsoo 



gsi^\r) = J2^2l + l)\Riico\r)\' + J^i^l + l)mco\r)\' . (25) 

1=0 1=0 

Since we are interested in the correction to the shift due to nonzero spacetime curvature, 
we will compare to the Lamb shift in flat spacetime Am, which is given, in the same 
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notation, by [13 



UJ — OOq CO + UJq 



To characterize the correction induced, we introduce a ratio between the two 



F r) 



A 



M 



(26) 



(27) 



Here, both and are formally divergent. However, the divergence can be dealt with 
by introducing a cut-off factor 18|, |28|] or resorting to a fully relativistic approach where 
no cut-off is present 29|, l30|. We choose the former in the present paper for simplicity. 

then the 



28 



In fact, if the cut-off is chosen as the electron mass as suggested by Bethe 
results from the two methods agree. As a result, we expect the ratio F{r) to be cut-off 
independent. To calculate the relative correction, F{r), at any given r, we need to know 
function gg(w|r) which depends on the radial parts, Ri{ijj\r) and Ri{uj\r), of the scalar field. 
Analytic solutions to the radial equation Eq. f|T4|) are however hard to find, so we now resort 
to a numerical approach [3l|, |32|. In this approach, one expands the radial functions Ri{uj\r) 

Ts/r) at a given position r outside the massive body 



as an infinite power series of (1 




1 ^iuir, 



where 



Suji{r) = ^afe(/, 



UJ 



k=0 



1 



r 



(28) 



(29) 



Inserting Eq. ( I28l) into Eq. (]T4l) . we obtain, after some simplifications, the recursive relation 
of the coefficients 



00 



— 2irs— j ak{l, cu) + 
ak-i{l,uj) + [3{k-2y 



CO 



-3{k - 1)2 + 2i{k - l)r 
+ /(/ + l)]afe_2(/, uj)-{k- 3) V_3(/, w) = 



c c 

2, 



(30) 



with ao{l,u) = 1 and afc(/,a;) = for /c < 0. For any given co and /, coefficients ak{l,uj) for 
any > can be found using the above relation. 

Now our task is to evaluate the transmission and reflection coefficients, Bi{u) and Ai{(jj), 
which can be determined by comparing the series representation of the radial functions fl25]) 
with their asymptotic forms (fT7|) and (fTSjl at large radii. For the radial function of the 
outgoing modes Ri{uj\r), one has, at infinity, i.e., r — )■ oo. 



After applying the operator 



^ d Lo r 
T> = - + — 



dr* 



(31) 



(32) 
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to both sides of the above equation, the right-hand side becomes 0(r ^). Then the reflection 
coefficient is found to approximately be 

Similar operations on the radial function of the ingoing modes Ri{u)\ r) yield the transmission 
coefficients 

2ir - 

- - KD-4k,W ■ ^''^ 

After simplifications, both coefficients can be expressed in terms of an infinite series with 
respect to ak{l,co) and (1 — r^/r) as 

Er=o[^(fc - 1) + - ^-4K(^, - 

and 

irs- 

In our numerical computation, these coefficients are evaluated for large and increasing r until 
they have converged to a desired accuracy such that |1 — + 15/(1^)1^)1 is less than a 

prescribed number Ai. Fig. [T] shows the numerically obtained values of |yl/(a;)p and \Bi{u)\'^ 
for various / and u for r = 5 x 10^ and Ai = 10~^. It is obvious that both coefficients 
change very rapidly with u. On one hand, our numerical result on the transmission coefficient 
is consistent with what one expects from a geometrical approximation, i.e., for a given 



quantum number I, \Bi{u)\'^ ~ 1 at high frequencies 25|], and on the other hand, it also 



agrees well with the approximate analytical expressions in the low frequency limit given by 



Page in Ref. [33 . 



Similarly, we can numerically evaluate S^i{r) for various u and / that range from zero to 
some large values such that a desired accuracy A2 is achieved. In this paper, we terminate the 



summation in S^i{r) for a given r, when 



< 10-1°. Now with S^i(r) and 



the reflection and transmission coefficients given, we can compute gg{uj\r), which is needed in 
our evaluation of A^, by appealing to Eqs. ( 128|1 and (l29l) . In the numerical computation of 
g^(ci;|r), we have set the accuracy to terminate the summation with respect to the quantum 
number / to be A^ = 10~^ such that — ^ ^ both the outgoing and 

ingoing modes. Fig. |2] shows how g^((X'|r) changes with u. Notice that g^(u;|r) is plotted 
instead of gg{u)\r) for convenience. One can see from Fig. [2] that g^(ci;|r) is a monotonously 
increasing function of u. 

The relative correction to the Lamb shift at different radii can now be obtained by 
numerically integrating both the denominator and the numerator in Eq. ( p7|) . Let us note 
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FIG. 1: and |i?;(a;)p as functions of u: for / = 0, 1,2,3,4,5. The horizontal axis is in the 

unit oi c/rg. The dashed hues represent |^;(a;)|^ and the sohd \Bi{ijj)\'^ . 
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FIG. 2: gg{ijj\r) as function of w, r = 5 r^, 10 r^, 15 r^. The horizontal axis is in the unit of 

c/Ts- 



that both the Lamb shift of an inertial atom in Minkowski spacetime Am and that of a static 
one outside a massive object are formally divergent. As we have already mentioned, we 
deal with the divergence by introducing a cut-off. For Am, the cut-off factor is chosen as 



where rrie is the mass of the electron [18|, |28 



However, for A^, a red-shift factor 
should be included. The reason is as follows. Our calculation of the Lamb shift involves the 
proper time of the atom (refer to Eqs. (jl]) and ^) but the field modes, Eqs. (IT^ and (IT^ . 
are expressed in terms of the Schwarzschild coordinate time. Therefore, the cut-off in A^ 

2 

should be chosen as ^/g^ since the frequency measured in the Schwarzschild coordinate 

'goQ. Taking the 



and that in the proper frame of the atom ujp are related by Ug 



energy- level spacing of the atom to be = 1.271 x lO^^s ^ which is comparable to the 



8 



transition frequency of a hydrogen atom and is within the frequency range of visible hght, 
we numerically calculate the relative correction to the Lamb shift of a two-level atom fixed at 
various radial distances in the unit of the Schwarzschild radius r^. We have computed cases 
corresponding to various masses of the massive body that range from 1.59 x 10^ ~ lO^'^kg 
and our result, which turns out to be independent of the mass of the massive body, is 
plotted in Fig. [31 This dependence should not come as a surprise since the radial distance is 
measured in terms of the Schwarzschild radius which varies with the mass of the object. 
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FIG. 3: The relative correction to the Lamb shift of a static atom at different radial position, r, 
which ranges from 1.1 r<j to 50 r^. 

A few comments are now in order. First, the relative correction F{r) is always smaller 
than unity. This means that at any position, the Lamb shift of the atom at a fixed radial 
distance outside a massive object is always smaller than that of an inertial one in flat 
spacetime. It is worth pointing out that this result seems incompatible with the analytical 
approximation ^ [1 + /(r)]AM previously obtained 34] (see Eq. (39) in Ref. js^) where 
/(r) is positive. This discrepancy is caused by the failure of estimating, in the asymptotic 
regions, the summation over the ingoing and outgoing modes, Ri{uj\r), Ri{(xj\r), to the same 
order of approximation in Eq. (29) and Eq. (30) in Ref. 3J] which are actually taken from 
Ref. j27|. In fact, to keep the approximation to the same order, terms which are of the 



order next to (x.^M/r) should be added both to Eq. (29) and Eq. (30) in Ref. 3J| (also 
to corresponding equations in Ref. 27|]). However, the analytical expressions for these 
terms are unfortunately hard to find. The numerical calculations we have performed here 
show that their contributions are actually negative and overtake /(r) so that the net effect 
is to make the Lamb shift smaller. As a result, /(r) in Eq. (33) in Ref. [34| should be 
replaced by a negative-valued function. We also want to point out that the Lamb shift 
we are calcula ting here is in close analogy with that near a conducting surface in a fiat 
spacetime [35|, l36j. In a fiat spacetime with a conducting surface, the vacuum field modes 
are completely reflected at the boundary, whereas, in the present case, they are partially 
reflected and this reflection occurs at every point in space in contrast to the flat spacetime 
case where the reflection happens only at the boundary. Second, near the horizon, the 
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Lamb shift of the atom decreases rapidly as the radial distance increases until it reaches 
the minimal value at ~ 2.0 where the correction is about 25%, then it grows very fast, 
but the slope flattens up at about 20 where the correction is still as large as 4.8%. After 
that the Lamb shift increases very slowly and is expected to approach the value in the 
Minkowski spacetime at infinity where the spacetime is asymptotically fiat. Remarkably, 
even for a radial distance as far as 50 r^, the correction is still as large as about 1.6%. 
Finally, let us note that r = (where the horizon is located) is a singular point, so in our 
numerical calculations, we can not start exactly at r = but rather at a point close to it 
and r = l.lr^ in Fig. [3] is just such a point. The closer this point is to the horizon, the 
closer the Lamb shift is to the value in the Minkowski spacetime. This is consistent with 
our previous result that the Lamb shift of the atom at the horizon approaches the value in 



the Minkowski spacetime [3^. As a matter of fact, the disappearance of the correction to 
the Lamb shift close to the horizon is a reflection of the fact that the effective potential, 
Eq. (fT6l) . which characterizes the scattering of the vacuum field modes off the spacetime 
curvature, diminishes to zero near the horizon (and at spatial infinity). This point is also 



discussed in Ref. 34 



So, for a very compact super-massive astrophysical body such as a neutron star or a 
super-massive black hole in the active galactic nuclei or even a stellar mass black hole, the 
curvature-induced correction to the Lamb shift is remarkably large. In principle, we can 
look at the spectra from a distant compact super-massive body to find such corrections. 
Therefore, our results suggest an interesting way to test fundamental quantum effects using 
astrophysical observations. In terms of the possible observation tests, one should keep in 
mind that our result is the Lamb shift as seen by a static observer at the position of the atom. 
For observation, the physically significant quantities are the ones related to the observations 
performed at a very large radial distance from the center of gravity. Let us note that the 
fact that the Lamb shift as observed by a distant observer will be gravitationally red-shifted 
makes the correction seem larger than it actually is. Table. [land Table. HTl show the relative 
correction to the Lamb shift near the horizon r ~ and at larger r respectively. Notice that 
F{r) and F'{r) represent, respectively, the relative Lamb shift observed by a static observer 
at the position of the atom and that by a distant one at the spatial infinity. 



TABLE I: The relative correction to the Lamb shift near r ~ r 



r/rs 


1.1 


1.3 


1.5 


1.7 


1.9 


2.1 


2.3 


2.5 


2.7 


2.9 


F(r) 


0.918 


0.821 


0.780 


0.759 


0.751 


0.752 


0.756 


0.759 


0.769 


0.777 


F'{r) 


0.277 


0.394 


0.450 


0.487 


0.517 


0.544 


0.568 


0.588 


0.610 


0.629 



Finally, let us point out that the correction to the Lamb shift we have just calculated is 
for an atom on a trajectory at a fixed radial distance outside a compact massive body. Now a 
question arises as to how the trajectory is realized, in other words, who is holding the atom, 
and this question is specially relevant when actual observational verifications are desired. A 
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TABLE II: The relative correction to the Lamb shift at larger r. 



rjrs 


5 


10 


15 


20 


25 


30 


35 


40 


45 


50 


F(r) 


0.840 


0.909 


0.939 


0.952 


0.964 


0.967 


0.976 


0.980 


0.983 


0.984 


F'{r) 


0.751 


0.862 


0.907 


0.928 


0.945 


0.951 


0.962 


0.968 


0.972 


0.974 



natural way to keep an atom at a fixed radial distance is to let it undergo circular geodesic 
motion. It has been shown that bound circular orbits are possible for massive particles only 
when r > AGM/c^ and among them only orbits with r > QGM/c? are stable |37|. The closer 
is the orbit to r = AGM/c^, the greater is the speed of the particle. However, an atom under 
geodesic circular motion will be subjected to the circular Unruh effect which has already 



been studied in the literature |38H41| and as such additional corrections to the Lamb shift 
will result. However, as we will demonstrate next, such corrections are in general much 
smaller than the curvature induced corrections for a typical very compact massive body 
such as a neutron star or a black hole. To this end, let us note that for a two-level atom in 
uniform circular motion, the correction to the Lamb shift is given in the high- velocity limit, 
v/c > 0.85, by jil 



Ac = ^^^^^ ^ (^e-2v^^-oc/aj.i^2 V35u;oc/a) - e^'^''^°'/''Ei{-2V3Buoc/a)^ , (37) 



where Ei denotes the principal value of the exponential integral function 42], a is the 



centripetal acceleration, uq the energy gap between the two levels of the atom, A = l + 

B = 1 — r(— )^ with 7 = (1 — ^)~^/^ and v is the velocity of the atom. Notice that here 



we have corrected a couple of typos in the original formula given in Ref. |41|. To estimate 
how large typically will be, let us take a neutron star of a solar mass M ~ lO^^kg as 
an example. If we choose r = 4.2GM/c^ which is very close to the inner most circular orbit 
and Uq ~ lO^^s"^ which is the typical hydrogen transition frequency, then we can easily 
calculate the velocity using the formula in Ref. 37| to get v ~ 0.913c and this leads to the 

2 2 

centripetal acceleration a = ^ 1.448 x 10^'^m/s^, A ^ 1.120 and B ^ 0.960. With these 
constants given, we can evaluate A^ and compare it with Am to find ^ ~ lO^^s. This 
ratio becomes even smaller when the radius of the circular orbit is larger. This shows that 
the corrections due to the circular motion to keep the atom at a fixed radial distance outside 
very compact massive objects is completely negligible as compared to those caused by the 
spacetime curvature. So, in practical sense, our results also apply to atoms orbiting a very 
compact massive astronomical object in circular motion at a given radial distance. 
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